Spin ice under pressure: symmetry enhancement and infinite order multicriticality 
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We study the low-temperature behaviour of spin ice when uniaxial pressure induces a tetragonal 
distortion. There is a phase transition between a Coulomb liquid and a fully magnetised phase. 
Unusually, it combines features of discontinuous and continuous transitions: the order parameter 
exhibits a jump, but this is accompanied by a divergent susceptibility and vanishing domain wall 
tension. All these aspects can be understood as a consequence of an emergent SU(2) symmetry at 
the critical point. We map out a possible experimental realisation. 

PACS numbers: 05.50.+q 05.70.Jk 75.10.Hk 75.40.Cx 75.40.Mg 



One fascinating aspect of condensed matter physics is 
the extent to which the nature and symmetries of emer- 
gent low-energy degrees of freedom can be independent 
of the high energy ones from which they derive. While 
many types of low energy degrees of freedom are possi- 
ble, systems in which these correspond to a gauge field 
are still rare. 

In this paper we study this theme in a setting which 
brings forth an unusual sequence of low energy fluctua- 
tions and symmetries: we consider spin ice, a frustrated 
Ising magnet on the pyrochlore lattice, whose low tem- 
perature behaviour is well-described by a gauge field, 
an emergent magnetostatics, representing a so-called 
Coulomb phase where correlations decay algebraically as 
r^^ . Transitions out of such a phase have attracted a 



great deal of attention recently 

We ask what happens when the host crystal is sub- 
jected to uniaxial pressure. This is a powerful probe of 
correlated matter, which has previously produced inter- 
esting information in frustrated magnets closely related 
to spin ice [HI, Q. Specifically, we consider the case of 
pressure-induced strain that lowers the crystal symme- 
try from cubic to tetragonal. We find that this induces a 
highly unusual symmetry-breaking phase transition out 
of the Coulomb phase: it is characterised by discontinu- 
ities in quantities such as the magnetisation, without be- 
ing first order. Rather, it occurs at an infinite order mul- 
ticritical point and is accompanied by a divergent suscep- 
tibility. We demonstrate both these features using Monte 
Carlo simulations. Further salient properties of the crit- 
ical point include the absence of a domain wall tension, 
along with spin correlations that vanish in a plane per- 
pendicular to the strain-induced tetragonal axis. These 
should be observable experimentally, and we discuss the 
scope for realising such a transition in the laboratory. 

The unusual nature of the transition is a consequence 
of an emergent SU(2) symmetry: as the pressure reduces 
the symmetry of the external space, an enhanced internal 
symmetry appears. We show this by means of an exact 



solution in three dimensions at the critical temperature 
and a mapping to a quantum phase transition in 2 -I- 1 
dimensions, between Ising and XY anisotropy in a spin- 
1/2 ferromagnet. 

Transitions with analogous features have been studied 
previously in the context of ferroelectrics 01, using ver- 
sions of the model for potassium dihydrogen phosphate 
(KDP) introduced by Slater Q, which in its two dimen- 
sional form is equivalent to the 6-vertex model and ex- 
actly solved 0. Spin ice has the potential to provide a 
much cleaner realisation of this physics than the transi- 
tion in ferroelectrics, since its magnetic degrees of free- 
dom are much more accurately represented by Ising vari- 
ables supported on a rigid lattice. 

Spin ice is well modelled by Ising spins parallel to 
their local easy-axis on the pyrochlore lattice [13] . At 
sufficiently low temperature we can limit ourselves to 
the highly degenerate ground state ensemble with two 
spins pointing in and two out of each tetrahedron (ice- 
rules) and consider effective nearest neighbour interac- 
tions [m [l3l . After coarse-graining the magnetisation, 
the discrete ice-rules can be written in terms of a con- 
tinuous divergence free condition from which emerge the 
algebraic correlations 13, 14 1 . Our goal is to lift the de- 
generacy of the resulting Coulomb phase via the exchange 
modulation 
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where Jij = J — S > for bonds on (001) planes per- 
pendicular to the strain axis, and Jij — J otherwise (see 
Fig. [T]). For (5 > 0, the Z2 symmetry of this Hamilto- 
nian is spontaneously broken at temperature Tc ^ 6 in 
favour of a state magnetised parallel to the [001] direc- 
tion. If the strain were very large, so that Tc ^ J, or- 
dering would be from a conventional paramagnetic state 
and in the standard three-dimensional Ising universality 
class. In contrast, for realistic, small strains is much 
less than J and ordering is from the Coulomb phase, with 
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the striking features we describe. 
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FIG. 1: A tetrahedron of the pyrochlore lattice with exchange 
J and J — S as indicated: the dashed bonds lie in (001) planes. 
One of the two ground states for 5 > is shown. 

Wc first present results from Monte Carlo simulations, 
done with a worm algorithm 15 , [l6| in which loop up- 
dates preserve the divergence free nature of the magneti- 
sation. The temperature dependence of the magnetisa- 
tion AI and of the inverse susceptibility are shown in 
Fig. [21 upper panel: remarkably, M takes its saturation 
value Msat at all temperatures below Tc, but x diverges 
as Tc is approached from above. To better characterise 
the transition, we also examine the probability distribu- 
tion function (PDF) of the order parameter, shown in 
Fig. [21 lower panel. This is a standard diagnostic: two 
peaks are expected in the PDF at the transition point 
when this is first order, but only one peak if it is contin- 
uous [l7| . Again, we find unconventional behaviour: a 
PDF that is uniform over all values of M at Tc- 

An exactly uniform order parameter distribution arises 
within Landau theory as a limiting case. The free energy 
density G{M) close to a multicritical point of order n is 



G(Af) = ^(T- 



2n' 



(2) 



In the limit n — > oo, it is independent of M at T ~ Tc- 
Moreover, the order parameter exponent takes the value 
/3 = l/(2n — 2), generating a jump for n — >■ oo. Also, the 
susceptibility {x A±\{T - Tc)/Tc\-^ for T ^ T^) has 
an exponent 7=1 independent of n and an amplitude 
ratio A_/A^ = l/(2n — 2) that vanishes as n — >■ 00. All 
these are borne out by our simulations. This success of 
mean field theory refiects the exceptional value of the 
upper critical dimension, dc2 = 1; which follows from 
standard arguments with allowance for the anisotropic 
spatial scaling, see Eq [5] below. 

In previous discussions of the ferroelectric KDP 
model 0, 0, 18-Ell| some distinctive aspects have been 
recognised. Indeed, Slater's original approximation Q 
yields a free energy independent of the order parame- 
ter at the transition, which on this basis has been called 
infinite order j2^ . Such a conclusion is evidently very 
delicate [31, and corrections to the approximation or 
physical perturbations have the potential to convert the 
transition into a conventional one: either first order or 
continuous. This is in fact the case for the material KDP 



itself, which has a first order transition that can be driven 
through a tricritical point under pressure Q. While in 
two dimensions, an exact solution of the KDP model cor- 
roborates the main features of Slater's results jol. [2ll. [23j . 
it is not a priori clear - not least in view of the famous 
difficulties in generalising Onsager's solution of the Ising 
model - what this tells us for higher dimension. Here we 
present simulations and exact results for the transition 
in three dimensions. In addition, we show how to realise 
it in a magnetic material. 
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FIG. 2: (Color online) Top: Simulation data for M (blue 
circles) and (x^c)"^ (red triangles) vs. T/Tc. Bottom: PDF 
of M/Msat at different temperatures. For T < Tc (red) the 
PDF has two sharp peaks at M = ±Msat- For T = Tc (black) 
the PDF is broad in small samples (dashed black line for 
L± = Lz — 4:) and almost flat in larger samples (solid black 
line for L± = A and = 18). For T > Tc (blue) the PDF is 
Gaussian. (Values at T ^ Tc scaled to flt on vertical axis.) 

The ice rules impose the same value of the magneti- 
sation M in all (001) planes: any configuration can be 
mapped onto an ensemble of strings of down spins span- 
ning the system from top to bottom, with the convention 
that the string vacuum is the state in which all spins are 
The number of strings defines (i) the total mag- 
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netisation and (ii) a topological sector, so that the former 
can be used to label the latter. The equiprobability of all 
sectors (see Fig. [2]) already suggests the absence of inter- 
actions between strings. This observation can be made 
precise using the transfer matrix that acts between two 
adjacent (001) layers of the lattice Since this trans- 
fer matrix is a direct product of factors T representing 
separate tetrahedra, we first discuss a single tetrahedron. 
A single string can enter the tetrahedron at either site in 
one (001) layer and leave it at either site in the other 
layer, with energy cost 4^/3. This sector is therefore rep- 
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resented in T by a 2 x 2 block in which all entries are 
the Boltzmann factor F = e"^''^/'^. By contrast, zero or 
two strings in a tetrahedron cost no energy and impose 
the next configuration locally, generating two 1x1 blocks 
with unit entries in T, which hence has the form 

r r 

T = p p . (3) 

V 1/ 

The maximum eigenvalues of T lie for T < Tc = -^j^ in 
the two fully magnetised sectors, but for T = T^. all sec- 
tors have the same maximal eigenvalue. In addition, the 
associated eigenvector in the one-string sector is (1,1). 
Explicit calculation of the transfer matrix for a complete 
lattice is complicated by the fact that the repeat unit 
in the [001] direction involves four different layers. Nev- 
ertheless, its maximal eigenvalues inherit precisely the 
properties described for T, and at T = Tc the asso- 
ciated eigenvectors in each sector give equal weight to 
all arrangements of strings on a (001) layer. The max- 
imal eigenvalues determine the physical properties for a 
system of size x L\ if the [001] dimension satisfies 
Lz S> L\. As a result: (i) below Tc the magnetisation is 
saturated; and (ii) at all sectors are equiprohahle and 
all configurations within a sector are equiyrobable. reflect- 
ing the absence of interactions between strings [1^1. By 
contrast, for a system with Lz ~ L± subleading eigen- 
values also contribute and the PDF for M is rounded, as 
is apparent in Fig. 

The transfer matrix for the complete lattice can alter- 
natively be thought of as the evolution operator in imag- 
inary time for a two-dimensional quantum XXZ model 
(in analogy with Eq. (8) of 0) with Hamiltonian 

'HQ = -jY^{st^ + sts^ + ^sts^ . (4) 

The strings thus denote the world lines of sf = 1/2 spins. 
In this language the decomposition into sectors is a con- 
sequence of the conservation of total s^, and is due to a 
U(l) symmetry. The ground state of H has ferromagnetic 
order of the quantum spins for all A, with an orientation 
for the magnetisation that depends on A. For A < 1, 
spins lie in the x-y plane and the state has zero-point fluc- 
tuations, representing the Coulomb phase. For A > 1, 
spins are aligned along the z-axis and the state has no 
fluctuations, representing the low temperature phase of 
the classical system. The equiprobability of string sec- 
tors at the critical point corresponds here to the degen- 
eracy of the N + 1 ground states of the isotropic Heisen- 
berg ferromagnet with TV spins, for all values of the total 
magnetisation = {—N/2, .., A^/2}. We have confirmed 
that this isotropy is a property of the full transfer matrix, 
and not only of the leading eigenvectors, by checking that 
the matrix commutes with the total spin raising and low- 
ering operators, = ±isp. This constitutes an 
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FIG. 3: (Color online) Mam panel: Behaviour of z ■ C{v,z) 
vs r'^/z, on a semi-log scale for 2 = (x), 2 (•), 4 ( ), 6 
(a), 8 (-(-), 10 (t). Dashed line shows behaviour expected 
from Eq. ((5]). The system size is L±_ — 14 and — 196. 
Inset: Same data plotted as C(r, z) vs r/L± on a hnear scale, 
showing that correlations vanish within (001) planes. 



enhancement of the symmetry at the critical point from 
U(l) to SU(2). 

The quantum description can be employed to calculate 
correlation functions within a given sector in spin ice, by 
approximating the imaginary time direction as contin- 
uous and treating Hq using harmonic spinwavc theory. 
For a system of shape sufhciently anisotropic (Lz ^ L\^) 
to justify a ground state treatment of the 2-1-1 dimen- 
sional problem, and taking z and r to denote distances 
along the [001] direction and in the (001) plane respec- 
tively, with p a microscopic length, two-point correlations 
at Tc are 

C(r,z)cx - expf-^^ . (5) 
z \ pz J 

As the strings do not interact, this form reflects the string 
auto-correlations present in a random walk in two dimen- 
sions with propagation time z. It agrees well with results 
from simulations, as we show in Fig. [3l 

A further consequence of the absence of interactions 
between strings is that the surface tension between op- 
positely magnetised domains vanishes and the domain 
wall width diverges, as Tc is approached from below. 
We investigate this phenomenon in simulations by us- 
ing weak, position-dependent magnetic fields to induce 
two domains with the interfaces between them lying on 
average in the x-z plane. Results are shown in Fig. |4l 
At Tc the magnetisation profile M{y) is a function only 
of y/L^_. Below Tc a mean field treatment of T-Lq gives 
a domain waU width ^"^ cx = (1 - T/TcY/"^: this is 
confirmed in Fig. [31 up to a constant vanishing in the 
thermodynamic limit. It is striking to find broad domain 
walls in an Ising magnet; they may be detectable using 
small angle neutron scattering. 

The natural way to generate the required degeneracy 
lifting in a magnetic compound is to apply uniaxial pres- 
sure along the [001] axis of a single crystal. We require 
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FIG. 4: Left: Magnetisation profile vs y/L± at T = 0.999 
for Lx = {6, 8, 10, 12, 16, 20}. Right: Domain wall width £^ 
vs \/t for L± = 20. The straight line is a guide to the eye. 
All samples have Lz = L\. 



6 positive. Within the nearest neighbour approximation 
the effective coupling has two contributions, one from the 
long-range dipolar interactions 11,12, 2^ , the other due 
to superexchange. While a uniaxial compression along 
[001] increases the former, the change of the latter is less 
clear as it depends on the evolution of orbital overlaps. 



Whereas this temperature is comfortably within reach 
of cryogenics, the dynamics of spin ice slows down greatly 
below Tf ^ 600mK, as activated ice-rule violating defects 
disappear. Their presence thus appears indispensable in 
practice, even though they will lead to more standard 
second-order phase transition in the three dimensional 
Ising universality class 0, [l^ ■ The width of its crit- 
ical region will, however, vanish with defect concentra- 
tion. Experiments should thus aim at the lowest dynam- 
ically accessible transition temperatures, with the target 
Tc 600mK only a factor of 3 away from existing ex- 
periments [i^ . We hope our work will stimulate further 
experimental efforts to realise this unusual transition. 
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under Grant PESC/RNP/HFM (PCWH and LJ), from 
ANR under Grant 05-BLAN-0105 (LJ), and from EPSRC 
(JTC) under Grant EP/D050952/1. 
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FIG. 5: Percentage change in magnetisation measured under 
uniaxial pressure of 10.5 kbar, M{P, B), compared to that at 
atmospheric pressure, M{0,B), at T = 1.7K as a function of 
applied field B along the [001] axis (extracted from [2^]'). 

Results from the one experiment so far performed on 
Dy2Ti2 07 under pressure [26| suggest that the effective 
coupling is indeed modified so as to produce S > 0. Pres- 
sure causes a 4% increase in M at small field, but a 1% 
decrease at high field (Fig. [5|) . The increase in the zero 
field susceptibility is consistent with the lifting of the de- 
generacy in favour of the states with magnetization along 
the field axis, as required for the transition studied here 
(Fig. [1]). The reduction of the saturated moment is ex- 
pected as the crystal is squeezed along [001] because of a 
resulting tilt of the easy axes away from [001]. 

In fitting the parameter S, we account for demagneti- 
sation effects by modelling the sample as a prolate el- 
lipsoid with major axes given by the dimensions of the 
approximately rectangular parallelepiped_specimen. For 
the interaction parameter of Dy2Ti207 
shown in Fig. [H yields - 200 mK. 



271, the best fit. 
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